Ferro-Orbitally Ordered Stripes in Systems with Alternating Orbital Order 
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We establish a novel mechanism of stripe formation in doped systems with alternating t2g orbital 
order — the stripe takes the form of a ferro-orbitally ordered domain wall separating domains with 
staggered order and allowing for deconfined motion of holes along the stripe. At a finite level of 
hole concentration this gives rise to the stability of this solitonic type of stripes, while we show that 
the phase change of the staggered order by tt plays a minor role in orbitally ordered systems. These 
results shed new light on the physics of doped materials in which orbital degeneracy is present. 
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In the context of superconductivity in cuprates both 
experimental and theoretical aspects of stripes have been 
the subjects of intensive research [Ij. At least in the 
low doping range, these systems can be viewed as an an- 
tiferromagnetic (AF) phase into which holes have been 
injected. When a hole hops in an AF background, it in- 
terchanges position with a single spin in each step and 
creates a "string" of flipped spins along its path ^] , which 
accumulates energy cost and thus confines the hole in the 
system with classical (Ising-type) interaction. At finite 
doping spin and charge density modulations (stripes) de- 
velop, which is a way to find compromise between two 
opposite tendencies: (i) to delocalize holes and gain hop- 
ping energy ex t, and (ii) let AF correlations to develop, 
which optimize the superexchange energy oc J. It seems 
that due to the presence of quantum spin fiuctuations 
stripes in cuprates show bond order, i.e. take the form of 
ladders with dominating singlet correlations on the rungs 
[3| . In addition, stripe formation has been demonstrated 
to exist in a model with classical AF exchange interaction 
[J|, so called t-J^ model, outlined in Figs. [Tfa) and[Ub). 
While the t-J^ model is hard to realize in spin systems, 
we show in this Letter that a related mechanism of stripe 
generation would work in t2g orbital systems. 

Some previous theoretical analyzes of stripe formation 
in orbitally degenerate systems lead to the conclusion 
that lattice distortions are essential for this kind of or- 
dering in systems with eg orbital degrees of freedom [5| . 
Stripes in pnictides with active i2g orbitals were sug- 
gested only very recently |6|] . We will demonstrate that in 
systems with t2g degeneracy and for large onsite Coulomb 
interaction U merely the interplay between hopping and 
the orbital superexchange interaction gives rise to the 
formation of ferro-orbitally (FO) ordered stripes as do- 
main walls (DWs) between regions with the alternating 
orbital (AO) order. With the goal of analyzing this phe- 
nomenon in detail we concentrate on the recently intro- 
duced strong-coupling version of the multi-orbital Hub- 



bard model for spinless fermions [t\ on the square lattice 
(when the spins form a ferromagnetic order). This model 
is applicable either to transition metal oxides with active 
t2g orbitals (when the tetragonal crystal field splits off 
the xy orbital from the {yz, zx} doublet filled by one 
electron at each site, as for instance in Sr2V04 3]), or to 
cold-atom systems [9\ with active p orbitals [IC | . 
The strong-correlation limit of the model reads: 
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FIG. 1: (color online) The mechanism of stripe formation 
in the spin Ising model: a hole doped at the DW kink (a) 
moves together with the kink (b). In contrast, a hole doped 
at the DW kink in the orbital t2g system (c) is confined to 
two sites. Two domains with AF (AO) order are shown by 
arrows (boxes); broken bonds are marked by x. 
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Here o and b refer to two t2g orbital flavors J7| , yz ^ a and 
zx = b, and the summations are carried over i d ab sites 
in the ab plane. The orbital superexchange J = At'^ /U 
and the effective next nearest neighbor hopping r = t'^ /U 
apply when C/ ^ t [ffl]. The pseudospin operator is Tf = 
5 {nia — noi), while the projection operator V removes 
from the Hilbert space states in which any site is doubly 
occupied; for more details see Ref. [7[. 

Single hole motion is in principle confined in the spin 
t-Jz model due to the potential well effect caused by 
the formation of strings. An effective way to avoid the 
string effect is to form an antiphase DW between two 
AF domains, consisting of two semilines separated by 
a transversal kink, and to create a hole at one of two 
sites nearest to the kink in that DW, see Figs. [Ha)-(b). 
The hole may be shifted along the wall without increas- 
ing the number of broken bonds. As the energy gain 
of —t is typically larger than the energy loss Jz/2 due 
to a broken bond, at a certain level of hole doping the 
energy decrease induced by free hopping of holes along 
the stripe will compensate the increase of the magnetic 
energy caused by the creation of the DW [4[ . 

Despite the similarity between the spin t-Jz model [J] 
and the t2g orbital model given by Eqs. (1-5), the mech- 
anism of stripe formation based on soliton-like motion of 
the kink-hole complex is not applicable to orbitally de- 
generate systems. This can be understood by analyzing 
the DW shown in Fig. [IJc). Boxes aligned along the a (b) 
direction represent b (a) orbital flavors in two domains 
with opposite phases of the AO order. Again, similarly 
to the Ising AF state, we have created an antiphase DW 
with a kink, and we have removed an electron from one 
of two sites nearest to the kink center. The downward 
shift of the hole by one lattice spacing is blocked by the 
b orbital below it. This follows directly from the form 
of the hopping ^, which is one-dimensional (ID) in or- 
dered ab planes 7]. The upward shift of the hole by one 
lattice spacing is allowed, but after that move the hole 
will be blocked again from above. Therefore, the hole 
and kink motion are confined in the orbitally degenerate 
system with the straight antiphase DW and a single kink 
in it [Fig. [ijc)]. Nevertheless, since the term ([2]) brings 
the biggest energy scale {t ^ J) we may expect that it 
will modify the form of the ground state above a certain 
filling level. The mechanism of hole deconfinement, how- 
ever, is different as we show below by a detailed stability 
analysis. 

The DW depicted in Fig. [l]is the most favorable one 
in terms of the minimal number of broken bonds per one 
DW site. For an orbitally degenerate system, however, 
the hole motion by hopping t is allowed along the chain 
only when orbitals reorient, similarly to the ID Cg sys- 
tems [12| . Such a FO ordered vertical chain of a orbitals 
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FIG. 2: (color online) Processes governing hole propagation 
in the stripe formed as a DW between two domains of the AO 
order (a): a doped hole may hop either along the DW (b)-(c) 
or sidewards (d)-(e); (f) shows the PW of a orbitals formed 
within a single AO domain; broken bonds are marked by x . 



in Fig. ^a) provides a DW between AO domains, and 
makes it possible to deconfine the hole motion along it. 
The price one has to pay are two (not one) broken bonds 
per site. After replacing one electron (orbital) in the 
chain by a hole [Fig. [2tb)], the hole can move due to Ht 
^ by one step [Fig. [D^c)] , and two broken bonds are re- 
moved and two other ones are created — hence the total 
number of broken bonds remains unchanged. In this way 
the hole motion occurs in both directions, leaving behind 
the undisturbed zig-zag pattern of broken bonds. 

The hole may also penetrate into the orbitally ordered 
domains, as for instance when the term ([2]) interchanges 
the hole with &-orbital in Fig. [UJb), leading to the state 
depicted in Fig. [2Jd). Now the hole can hop further to a 
nearest neighbor site either downwards or upwards, but 
when it enters "deeper" into the AO domain, the num- 
ber of broken bonds increases (now by one) and the string 
effect [4I occurs [Fig. [He)]. This mechanism efficiently 
confines the hole motion to the stripe DW. The expected 
energy gain due to hole motion over a homogenous AO 
state is realized by the parallel alignment of a-orbitals 
in the chain, while the phase change by tt between two 
AO domains is irrelevant to achieve that gain. A similar 
FO ordered chain deconfining the hole motion and hav- 
ing the same energy cost can also be created by reversing 
every second orbital along a vertical polaronic line within 
a single domain of the AO order, see Fig. ^i) — it gen- 
erates a polaronic wall (PW) by a mechanism similar to 
that which operates in the ID eg orbital model [l2J. The 
main difference between stripes depicted in Figs. [2Kb) 
and[2Kf) is that the hole in the DW stripe may penetrate 
one domain from each stripe site, while in the PW stripe 
it can enter both domains from every second site [the hole 
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FIG. 3: (color online) Artist's view of hole propagation in the 
homogenous |AO) state (a) — a hole may move due to t @ 
and create broken bonds, marked by x, after one (b), two (c), 
or three (d) hops, while r (|4]) does not disturb the AO order 
(a) and (e), but may also leave behind broken bonds (c) and 

(f)- 



in Fig. ^i) can move only from site j (not from I) cither 
to site i or to site k]. The explicit analysis presented 
below will demonstrate that it is energetically somewhat 
more favorable to create the antiphase DW stripe. 

In order to discuss stripe stability we must first analyze 
single hole propagation in the homogenous AO phase by 
means of the same method which will be used later for 
stripes. We begin with a hole replacing a orbital [Fig. 
[3lja)] — the hopping term oc t can shift a hole only hor- 
izontally to its neighboring sites which are occupied by 
b orbitals [Fig. [31[b)]. The hole motion can be contin- 
ued and after the second and the third hop the number 
of broken bonds increases, see Figs. |3fc)-(d). In gen- 
eral, the energy rise oc J, growing with the number of 
broken bonds, confines the hole motion (here initiated 
by the hop in the direction a) [3|- The number of de- 
fected bonds can be easily evaluated explicitly for paths 
up to length of seven steps, as we have done in our anal- 
ysis. For longer paths reasonable approximations can be 
developed. They follow from the relation between the 
number of broken bonds and the number of bends in the 
path, see Fig. ^d), and on the number of zigzags, such 
as that arising if the hole in Fig. ^c) has moved left. 

The action of the term ^ does not give rise to hole 
deconfinement, as it moves a hole generating several bro- 
ken bonds, cf. Figs. [3l^a)-(c). In contrast, weak hole 
deconfinement (r <^ t) occurs due to the term ^ which 
shifts a hole by two lattice spacings. This takes place 
without bringing about any additional defects in the AO 
order and is allowed provided that the opposite orbital 
(here b) occupies the intermediate site [Figs. [31Ja),[31Je)]. 

We are going now to cast the insights which have 
been outlined above into the framework of the recursion 



method 13| applied to the Green's function, and to de- 
termine the self-energy. The starting point is the bare 
Green's function Go(fcb,w) related with hole deconfined 
movement in the |A0) state mediated by the free prop- 
agation term (J4l) along the b axis, with momentum kb, 
[Fig. Etc)] — it is Go^(kb,u;) = cu ~ 2TCOs(2fc&) - J. A 
reference energy of the |A0) state has been subtracted 
in GQ^{kb,uj), and the energy J above arises from the 
four bonds removed from the jAO) state by adding a 
hole. The energy dispersion oc r in Go(fcf,,w) is given by 
the matrix element of the full Hamiltonian evaluated for 



the propagating state |1) = 



X;„exp(i2nfc, 
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where L is the system length along b, and the hole is at 
the origin in Fig. ^a). 

The full Green's function contains the selfenergy, 
G^^{kb,u}) = GQ^{kb,uj) — S(a;), which stems from the 
confined motion initiated by the first step in the a direc- 
tion transverse to the coherent propagation along b. The 
confined motion [Figs. [3l[a)-(d)] is accompanied by string 
formation and path retracement by holes. We evaluated 
I](a;) analytically by applying the recursion procedure, 
i.e. by the consecutive action with the Hamiltonian ([T]) 
on the state [1), which represents the hole created in the 
perfect AO order [Fig. [SJa)]. The only approximations 
consist of neglecting: (i) all details in long hole paths, (ii) 
some processes mediated by the term ^ (such as those 
represented by Figs, ^c) and[3ff)rthey bring about only 
a minor incoherent contribution [7|), and (iii) hole-hole 
interactions within AO domains. On the other hand, we 
always implement the constraint that the hole path can 
not encircle a plaquette — its topological reason can be 
recognized by analyzing Fig. [3l 

'S{uj) takes the form of a continued fraction. For ex- 
ample, if we neglect, for demonstration purposes only the 
term ([5]) (it is taken into account in the actual calculation 
but gives only minor corrections as r <C i), one finds 
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where prefactors 2 at t^ refer to the number of ways 
by which a path can be further extended after the 1st 
and 2nd hop to a neighboring site. Such terms represent 
off-diagonal matrix elements of the Hamiltonian between 
consecutive states created during the recurrence proce- 
dure 13|. The prefactors 7 and 9 at J/4 stand for the 
number of bonds for which the superexchange ([3]) gives 
for paths of length 1 and 2 respectively (instead of —J/4 
as for all bonds in the |A0) state) — they enter as diag- 
onal matrix elements of the Hamiltonian. By searching 
for zeros of G^^ (fcf,, uj) the quasiparticle dispersion eo{kb) 
can be derived. For the numerical results shown in Fig. 
HJa) sufficiently long paths (up to 12 steps) have been 
considered in order to obtain the saturation of results. 

The same method of analysis apart from a minor 
modification may be also applied to the DW stripe it- 




phases using the energy gain per site {X = DW,PW), 



FIG. 4: (color online) (a) Quasiparticle dispersion tx{kb) for 
the homogeneous AO phase (X=0, dotted) and for two kinds 
of stripes: X=PW (dashed) and X=DW (solid line), (b) 
Energy gain SEx{ni) per site ^ for two stripe phases for 
increasing stripe filling nj. Parameters: J — OAt, r = O.li. 



self. Since two states shown in Figs. [2lb) and [Ifc) 
are not related to each other by a translation, both of 

them, |1') = A/|E„exp(z2nfcb)a2„^|DW) and |1") ^ 



rE„exp{*(2n+l)fcb}c 



'(2n+l)b 



|DW), need to be cho- 



sen as a starting point in the recursion method. Here 
DW) stands for the empty stripe depicted in Fig. ^a.) 
and the position of the hole in Fig. EJb) is given by 
the zero vector. Thus the recursion method has to be 
generalized to the case of several initial states. We have 
found this generalization, but since the procedure has 
turned out to be e quiv alent to the so called projection- 
operator technique [l4|, we do not discuss it here in de- 
tail. The Green's function and the self-energy are given 
now by 2 X 2 matrices. The inverse bare Green's func- 
tion has only off-diagonal nonvanishing matrix elements, 
[GQ\h,io)]vv, = [G(7^(fcb,w)]i"i' =Lo-2tcos{h)~J/4:, 
where the reference energy due to the Ising term ([3|) has 
been subtracted again. When applying the recurrence 
procedure to the evaluation of S(a;) we may use our pre- 
vious observations regarding the matrix elements of the 
Hamiltonian between states representing different paths, 
in some cases including modifications brought about by 
the presence of the DW. 

The PW stripe formed as a FO ordered chain within 
the AO order, shown in Fig. [Hf), can be also solved 
by using the generalized version of the recursion method. 
By looking for the zeros of the inverted Green's function, 
or of the determinant, if G{kh,io) is given by a matrix, 
we have determined and showed in Fig. lU^a) the energy 
dispersion £o{kb) in the b direction for the homogeneous 
system with AO order, epw(fcb) for the PW stripe, and 
ejyw{kii) for the antiphase DW stripe. Different disper- 
sions confirm that the hole motion stems from t (t) hop- 
ping in the stripe phases (homogeneous AO phase). 

We investigate the stability of both types of stripe 
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dkbex{kb)-Eo{ni) + 5Ej, (7) 



with respect to the doped homogeneous AO phase with 
energy Eo{ni) — n;eo(7r/2), given by the band energy 
minimum eo(7r/2) and proportional to the linear stripe 
filling ni in the low doping regime (when the linear filling 
does not increase the global filling). Here SEj = J/2 is 
the energy cost of two broken bonds in a stripe phase. 
The stripes are stabilized by increasing ni, see Fig. IDJb), 
but the PW stripes are somewhat less stable which we 
interpret as following from the destructive interference of 
hole penetration paths into left/right AO domains. 

In conclusion, we have shown that a purely electronic 
mechanism leads to self-organization in the form of FO 
ordered stripes at antiphase DWs, penetrating into the 
AO order — this novel phase becomes more stable than 
the doped homogeneous AO state at the linear filling of 
ni ~ 0.26 {ni ~ 0.25) for J = OAt (J = 0.2i). The energy 
gain over the PW stripes is small for all ni , which suggests 
that also the latter form of stripes might be formed at 
finite temperature. These features are unique and can 
be of relevance to the behavior of doped Mott insulators 
with t2g or p orbital order when spins may be neglected. 
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